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Lattice perturbation theory for gluonic and fermionic actions. 
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We calculate the two loop Landau mean links and the one loop renormalisation of the anisotropy for Wilson and 
improved SU(3) gauge actions, using twisted boundary conditions as a gauge invariant infrared regulator. We 
show these accurately describe simulated results, and outline a method for generating Feynman rules for general 
lattice field theories, in a form suitable for efficient numerical calculation of perturbative loop diagrams. 



The conflicting demands of controlling finite 
volume effects and regulating the computational 
effort are balanced by using a large spatial lat- 
tice spacing, Og. Symanzik improvement can be 
employed to control discretisation effects. The 
temporal lattice spacing, at, should remain small 
to give good resolution of correlation functions 
and control the finite at effects; improvement in 
the temporal direction would lead to 'ghost' poles 
in the gluon propagator. Whilst general direc- 
tional asymmetries are best couched in terms of 
a metric structure, we restrict ourselves here to 
the case where only the temporal lattice spacing 
is reduced. 

The anisotropy, x = '^s/cLt, is renormalised. 
Physical scales are related by the measured value, 
Xr = Z{g'^,x)Xi which we calculate in pertur- 
bation theory (PT). With tadpole improvement 
(using two loop (2L ) factors derived in Sec. 0), 
the one loop (IL ) determination of Z in Sec. ||is 
sufficient. 

1. The perturbative action 

We follow the notation of Liischer and Weisz 
|l) . The links of the cubical lattice in D dimen- 
sions are labeled (x, ^), ^ = 1, D and the basis 
of the unit cell is {e^}. 

The gauge field denoted U is J7,i(v) = 
exp{gA^{'v + where g is the bare coupling 

constant. A general action may be written as 
a sum of Schwinger lines made up of links. In 
the gluonic case, these are closed contours. For a 
fermionic action, the path is terminated at points 
X, y by quark fields ip, ip respectively. The per- 



turbative action and associated vertex functions 
conic from expanding the action as a polynomial 
in A, e.g. in the gluonic case: 

r ' Bi..Br 

with Bi = (ki, ^i, ai). The Feynman rule for 
the r-point gluon vertex is — V^. The con- 
cept is to use two separate programs; one is 
used once to produce data files encoding Vr 
from the action expansion. The other then 
uses these data files many times in evaluating 
loop integrals. Rather than store a very large 
lookup table for all possible colour, Lorentz and 
momentum structure we focus on the action- 
specific features. Writing Vr{Bi, Br) = 
Cr{ai, Gr) yr(ki, /ii; ...; kr, /Xr), the action 
independent colour factor is Cr(ai, a^) = 
Tr(T,,...r,J + (-l)'-Tr(T,„...r,J for a gluonic 
action. We include it only in the loop integration 
code. In addition we can write for a given action 
1^ (/x^ ; ki , .. ., k^) = 

"E'^-pfz;^k..vr(/.j). (1.1) 

Q=l ■ \ i=l / 

For given Lorentz indices, /i.^. = (^i, /i^), we 
have rir terms (clearly action dependent). Each 
term arises from picking a sequence of r gauge 
fields from the action with specific Lorentz in- 
dices, and the positions of each gauge field are 
the v's. For simple actions, nr{fjLr) = for many 
/Xj,. So, for given /x^, we can reconstruct Vr for 
general momenta and colour structure knowing 
only ({./■«;<, •■•,<} I a = 1, ...,nr)- 
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Figure 1. Feynman diagrams for the Landau mean 
link. E refers to the diagrams in Fig. ^. 



This must, however, be done r! times, as we 
must consider all permutations of the vertex legs 
in the loop integration code. We may offset this 
cost by carrying out some of the symmetrisation 
at the action expansion stage. This will be at 
the expense of increasing the nj.{pr) by the ad- 
dition of permuted terms. We do not consider 
colour factors in the action expansion, but a gen- 
eral permutation, ct, will change Cr- We thus 
break the symmetrisation down into two stages; 
the first comprises a that alter Cr by a factor in- 
dependent of the colour index values. These are 
the cyclic permutations and the inversion, where 
a-Cr = Xr(o')Cr, and can be included in the Tay- 
lor expansion code by absorbing XrCc) into the fa 
of the permuted terms. For r > 4, the remaining 
permutations are colour dependent and must be 
carried out in the loop integration code. 

We perform the Taylor expansion using a 
Python code, chosen for its object orientation 
and list handling properties. For a more detailed 
specification of the implementation see The 
Taylor expansion of the links in a path (or paths) 
is represented by an instance of the taylor class, 
F. This is implemented as a 'dictionary' whose 
keys are lists of Lorentz indices, fi = [fii, ...,/ir] 
for r up to some chosen maximum order. F[fj,] 
returns a list of rir instances of the entity class 
of order r. Each encodes a monomial contribu- 



tion to Eqn. (1.1) and may be implemented as a 
list: Er = [/, X, y, vi, Vr], where x, y are re- 
dundant for gluonic actions. To illustrate, for the 
single link U^{z): 



F[[]] = [l,z,z + e^] 
FM] = [l,z,z + e^,; 
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Figure 2. Feynman diagrams for the self energy 
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etc. It is useful to store position vector com- 
ponents as (integer) multiples of a/2. We de- 
mand that the product of two taylors represent- 
ing two paths should be the Taylor expansion of 
the combined contour obtained by following the 
first path, then the second. Gauge covariance 
demands that the second taylor be thus trans- 
lated to start where the first path ceased. In the 
above example, F' = F * F \s the expansion of 
Uf^{z)Uf^{z + e^), and F' * F' would be a four 
link object. This allows for compact definitions of 
actions using templates that are themselves com- 
posite, or recursive, functions of other objects. In 
particular, covariant derivatives and fat links can 
be described without explicitly expanding all the 
contributing paths. 

For a pure gauge action, as well as being trans- 
lationally invariant, the final paths are closed and 
traced, and with only the real part taken. This 
high degree of symmetry permits identification 
of many equivalent terms in F, and a signifi- 
cant compression of the taylor can be carried 
out prior to symmetrisation. 

1.1. the propagator 

The gluon propagator is computed for given 
momentum by inverting the two-point function 

(2) (2) 

VJiV having fixed the gauge by adding Vi2 u an 
extra term. The inverse propagator is then 
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r(";(k) = v;(^Hk) 
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Figure 3. 2L Landau mean link coefficients for the 
Wilson action. 

where fc^ = 2 sin(fc^/2) and a = 1 is Feynman 
gauge. Inverting this, a general gauge 7 (such as 
Landau gauge, 7 = 0) has a propagator (up to a 
colour factor) 



G(x)(k)^G(";(k) 



ix 



21,2 



We follow § and use twisted periodic bound- 
ary conditions for the gauge field. There is no zero 
mode and no concomitant infrared divergences 
whilst gauge invariance is still maintained. 

The twisted boundary condition for gauge 
fields is J7^(x + L^e^) = 0i,[/^(x)0~\ where 
the twist matrices 0i, are constant SU(N) matri- 
ces which satisfy 0^0,^ — Zf^^O^Qf^ with 



exp(27rmpy/A^) an element of the centre of SU(N) 
with 71^1^ — 0, l,...,iV — 1, an antisymmetric 
tensor. In the calculations presented below we 
use 'doubly twisted' boundary conditions: ni2 — 

-^21 = 1, n 



otherwise, and also 'quadruply 
twisted': UijCijk = (1,1,1), noj = (1,1,1). The 
former is used as in ||] to define a physical gluon 
mode and the latter most effectively suppresses 
tunnelling between equivalent vacua and so allows 
the simulation of the Landau meanlink to match 
PT. It is important that nijnok^ijk = mod3 so 
that the instanton number is zero 

For the double twist with lattice of extent L in 
the 1,2 directions and continous in 3,4 the mo- 
mentum spectrum is continuous in k^, /c4 and dis- 
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Figure 4. 2L Landau mean link coefficients for the 
SI action. 
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i = 1,2, and ki,ni G Z. The interstitial mo- 
menta, rii, play the role of the colour index and 
excluding ni = 712 = eliminates the zero mode 
and imposes a gauge- invariant infrared cutoff mo- 
mentum of 2tt/NL. The Clebsch-Gordan factor 
is now a complex phase |^ . When a graph is built 
it must have the same overall centre phase factor 
as the tree diagram of the same external topol- 
ogy. For example, the self energy must have the 
same factor as the inverse propagator. 

The Feynman rules arising from the ghost fields 
used in the Fadeev-Popov gauge fixing, and from 
the measure are given in [Q . 

We consider the Wilson action (W) and the 
Symanzik improved action (SI) in Q, 

Sw{f3,x)^-P (xPs,t + -Ps.. 



Ssi{P,x,v)^~P[x\7;Ps.t 



12 V 



-P I — 



1 

12 



summing over spatial directions s, s' < s. P 
are plaquettes, and R denotes 2x1 loops. We 
use these forms of the actions as they minimise 
the dependence of the action parameters on the 
coupling, ~ &/ f3. 



Table 1 

Fits to perturbative expansion coefficients as functions of x, as defined in the text. 



action 


quantity 


constant 


1/x 




log(x)/x 


Wilson 




-0.01232 (2) 


0.0931 (1) 


-0.048 (1) 


-0.049 (1) 




t 




0.09132 (2) 


0.01385 (9) 


0.0090 (1) 




Ujs 


-0.017 (1) 


019 


-0 023 (4) 


-0.041 (4) 








0.0036 (2) 


-0.0246 (2) 


-0.011 (1) 




bi'^ 


0.030 (1) 


-0.042 (3) 


0.015 (4) 


-0.013 (4) 




b?' 




0.0103 (2) 


-0.0072 (2) 


-0.038 (1) 




V 


0.1687 (2) 


-0.16397 (4) 


-0.005245 (2) 




Symanzik 




-0.1012 (2) 


0.0895 (4) 


-0.0513 (6) 


-0.0502 (7) 


improved 


b^^ 




-0.0194 (1) 


-0.0453 (1) 


0.0110 (1) 






0.0091017 


-0.071095 


0.065993 


0.051785 




b?^ 


-0.007214 


-0.016691 


0.028623 


0.028351 




V 


0.0955 (4) 


-0.0702 (16) 


-0.0399 (14) 








0.0602 (1) 


-0.0656 (2) 


-0.0237 (1) 





Tadpole improvement factors, u^.t, arising from 
chosen self-consistency conditions amount to a 
reparametrisation of Sw 

In the SI case, we must also impose v — 
and PT is specific to a particular improvement 
scheme. If Ms = 1 + dsg^ + 0{g'^), then 

Ssi{P, X, v) = SsiiP, x,v = l)+ g^ASsi 

ASsi= ^ {xRs.,t + x~HRssy + Rs's'^s)) 

and we treat g^ASsi as a counterterm insertion 
(Fig. 1(f)) in the gluon propagator. 

2. Landau mean link improvement 

We focus on the Landau mean link scheme, 
where Us^t are defined as the expectation values 
of spatial and temporal links in Landau gauge. 

With6« =a«,5f) =aP)+a«(3ai^)+a«), 
the series for ui^s.t can be written as 

\l + b['^go'+bPgo'+Oigo'). 

The relevant diagrams are shown in Fig.^ The 
2L contributions were computed using a parallel 
version of Vegas. Fits to the L — > oo IL and 2L 



contributions are given in Table |l| and shown in 
Figs. I I 

By comparing such truncated series with re- 
sults from high-/3 Monte Carlo simulations, higher 
order coefficients in the perturbative expan- 
sion may be inferred Field configurations 
were generated using a 2nd order Runge-Kutta 
Langevin updating routine. The implementa- 
tion of the Langevin evolution is such that any 
pure gauge action can be simulated by simply 
specifying the list of paths and the associated 
couplings. The group derivative for each loop 
is then computed automatically by an iterative 
procedure which moves around the loop link by 
link constructing the appropriate traceless anti- 
Hermitian contribution to the Langevin veloc- 
ity. This is the most efficient implementation, 
minimising the number of group multiplications 
needed. It applies whenever the quantity to be 
differentiated is specified as a Wilson path. 

The twisted boundary conditions are imple- 
mented in the manner suggested in 0] where the 
field simulated, U = U everywhere, save J7^(x)0p 
when = L^, where 0^ is the twist matrix as- 
sociated with the fi direction and is the lat- 
tice extent. The action S{U) is identical to the 
untwisted action except that a loop whose pro- 
jection onto the (/^, t/)-plane, ^ < v, encircles the 
point (Lp + ^,Liy + i) has an additional factor 
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Figure 5. Deviation of the Landau mean link from 
2L PT as a function of go^ for the Wilson action. 

of (z^jy)"^ where c is the integer winding number. 
The program assigns the correct phase factors to 
each path appropriately. A Fourier accelerated 
algorithm is vital when fixing to Landau gauge, 
which is otherwise prohibitively time consuming. 

Both the mean link and plaquette values are 
measured for a range of (3 and x values. In Fig. |^ 
we plot for the isotropic Wilson action the sim- 
ulation results (with 2L PT already subtracted) 
versus g^^, showing the deviation from the 2L PT. 
Care must be taken when fitting polynomials to 
data in a given range of the expansion parameter. 
If the true coefficients are such that in this finite 
window there is an approximate cancellation be- 
tween some combination of terms, then all such 
terms must be included in the fit function. If not, 
as might be the case when trying polynomials of 
increasing order, marked instability in goodness 
of fit and even low order coefficients will be seen. 
Whilst this did not seem to be the case here, we fit 
the data points using 'Baysian' techniques B to 



obtain h 



(3) 



-0.00032 (5), h. 



(4) 



We believe that the value of bs is reliable and 
is plausible. 



-0.00027 (5). 



3. The renormalised anisotropy 

We calculate Z{g^ , x) = 1 + vix)9^ to IL from 
the dispersion relation for the gluon propagator 
in twisted QCD 

The calculations were done on lattices with 
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Figure 6. The IL anisotropy renormalization coeffi- 
cient versus anisotropy for the Symanzik action. 

SL^P twisted momentum grid points, where 4 < 
L < 32 is the lattice size in the twisted 1,2- 
directions and / > 50 gives an effectively con- 
tinuous momentum spectrum in the untwisted 
0, 3-directions. The gluon mode has momen- 
tum parametrised by p = {ipo,m,0,p3) , m = 
2tt/NL, where m is the pole mass of the bare 
gluon propagator. We calculate the pole energy 
of the propagator as a function of p^. For suf- 
ficiently small P3 the dispersion relation can be 
fitted to the standard quadratic form using the 
renormalised mass, m^', and XR parameters. 
The dispersion relation is derived from 



ri^,i";(p) = r(:)(p) 



To 0{g^) the momentum argument of I]^y(p) is 
given by the zero of the bare inverse propagator, 
most conveniently done for fi — 2. So, for given L, 
P3, we determine the bare pole value po = Eq by 
numerical solution of r2^''(p) = 0, and the bare 
mass mo is defined by Eq = mo/x at = 0. 

The lattice energy function is determined by 
the inverse propagator, F{E) = 4sinh^(_B/2). 
We define XR using dispersion relations for very 
small p3 

Fijno/x) , 
- F{mR/xR) , 



F{Eo) - pI/x'- 
FiEn) = pI/xI 
with Er defined by 



xF{Eo) ~ xF{En) = g^^{Eo, m, 0,^3) 



P(X=4) 




Figure 7. Comparing Z from simulation 
PT for the Wilson action. 



and IL 



To 0{g^) we define mn = mo + g^mi and 
V 



2x 



2x 



/^(L) and mi{L) are determined by a straight line 
fit, and extrapolated L oo. 

The IL diagrams for the self-energy Ep,y(p) are 
shown in Fig. ^j. For on-shell external momentum 
{iEo,m,0,p3), graphs Fig. ||(a,c) exhibit a pole 
in the integrand so the kg integration contour is 
shifted ko ^ ko — iEo/2. We also use the change 
of variables k ^ k' and fc^t = k'^ — a^sin(fc^). 
A good choice is 1 — (x^^)^^- The integrals 

did not require Vegas and took between 2 and 
16 hours on a PC. Feynman gauge was used. 

For the Wilson action we verify 77 (x — 1) = 
1 y L. For L 00, the mass renormalisation 
nil /mo agreed with and r]{x > 1) with pT| , 
both derived using the background field method. 
In Fig. we compare the IL PT prediction for 
Z with measurements from 17|] . The 2L mean 
link expressions from Table ^phave been used to 
rescale the coupling. The error in xr is < 5% for 
(3 > 5.8 and < 2% for /3 > 6.3 for 1.5 < x < 6. 
This is partly due to the error in the expressions 
for Us^t, e.g. Fig. ^, and the truncation error in 
Xr is expected to be smaller. 

Using the same method, for the SI action we 
consider actions with mean link improvement 



Table 2 

Comparing XR from simulation 
the SI action. 



and IL PT for 



/So 


Xo 


ajVMeV 


(mcas.) 
Xr 


Xr 


1.7 


4 


661 (11) 


3.56 (2) 


3.601 




6 


779 (28) 


5.28 (2) 


5.463 


1.8 


4 


797 (21) 


3.61 (2) 


3.629 




6 


839 (9) 


5.31 (2) 


5.489 


2.4 


3 


1200 (50) 


2.71 (3) 


2.776 



(ds = bi^^) and without (dg = and omitting 
Fig. 11(f)), obtaining 77 and rj' respectively. These 
are shown in Fig. ^, and fits to x are given in Ta- 
ble |. In Table | we compare these predictions to 
measured data from |^J^. The gauge couplings 
simulated are relatively small. For xo < 4 the 
error on XR is 2%, and for xo < 6 only 3%. As 
the simulated values of Us.t were used to rescale 
Xoi Po, these errors reflect the IL truncation error 
on Z(5o^x)• 

IL calculations of Z thus seem to be sufficiently 
accurate for use in setting the scale etc. in most 
lattice applications. 

We are pleased to acknowledge the use of the 
Hitachi SR2201 at the University of Tokyo Com- 
puting Centre and the Cambridge-Cranfield High 
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